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We develop the formalism to investigate the relation between the evolution of the large-scale 
(quasi) linear structure and that of the small-scale nonlinear structure in Newtonian cosmology 
within the Lagrangian framework. In doing so, we first derive the standard Friedmann expansion 
law using the averaging procedure over the present horizon scale. Then the large-scale (quasi) linear 
flow is defined by averaging the full trajectory field over a large-scale domain, but much smaller 
than the horizon scale. The rest of the full trajectory field is supposed to describe small-scale 
, nonlinear dynamics. We obtain the evolution equations for the large-scale and small-scale parts of 

■ the trajectory field. These are coupled to each other in most general situations. 

^\ ' It is shown that if the shear deformation of fluid elements is ignored in the averaged large-scale 

dynamics, the small-scale dynamics is described by Newtonian dynamics in an effective Friedmann- 
Robertson- Walker (FRW) background with a local scale factor. The local scale factor is defined by 
the sum of the global scale factor and the expansion deformation of the averaged large-scale displace- 
ment field. This means that the evolution of small-scale fluctuations is influenced by the surrounding 
qq ' large-scale structure through the modification of FRW scale factor. The effect might play an impor- 

tant role in the structure formation scenario. Furthermore, it is argued that the so-called optimized 
or truncated Lagrangian perturbation theory is a good approximation in investigating the large-scale 
^ ' structure formation up to the quasi nonlinear regime, even when the small-scale fluctuations are in 

OS ! the non-linear regime. 
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1. INTRODUCTION 



The observation of the anisotropy of the cosmic microwave background (CMB) indicates that the universe is 
remarkably isotropic on the present horizon scale. Thus it is natural to describe the horizon scale spatial geometry 
of the universe by a homogeneous and isotropic metric, namely, the Friedmann-Robertson- Walker (FRW) model. 
However, the real universe is neither isotropic nor homogeneous on local scales and has a hierarchical structure such 
as galaxies, clusters of galaxies, superclusters of galaxies and so on. It has been naively regarded that the FRW model 
is a large scale average of a locally inhomogeneous real universe. There have been several studies in this direction in 



general relativity |lj-|l6| MEMLI I 



Aside from such a fundamental problem, there is an interesting and practical problem associated with inhomo- 
geneities across various scales. Can one ask if the formation of small-scale structures is influenced by the gravitational 
effect of structures with larger scales? Such an environmental effect may be important and even essential to clarify 
the process of the hierarchical structure formation. 

This is the problem we attack in the present paper. Namely, we develop the formalism to investigate the gravitational 
instability in general situations where the large-scale linear and the small-scale non- linear fluctuations coexist. If one 
uses the N-body simulation to answer the above question, one needs high spatial resolution over a very large box 
comparable with the horizon scale, which may be well above the ability of the present computer. However, it seems 
reasonable to regard the situation such as local nonlinear structures are superimposed on a smoothed large-scale 
linear structure and the large-scale dynamics may well be treated by Zel'dovich-type approximations for an usual 
power spectrum. This suggests us to adopt an analytical approach based on the Lagrangian perturbation theory in 
Newtonian cosmology. 
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The reason why we consider Newtonian cosmology is partly because of its simplicity and partly because the New- 
tonian cosmology is a good approximation to a realistic inhomogeneous universe. In fact we have shown that the 
Newtonian cosmology in the relativistic framework is a good approximation even for the perturbations not only inside 
but also beyond the present horizon scale |[o],[5l) . 

The reason why we work within the Lagrangian framework is because it seems easy to introduce the averaging 
process. This is essential in our formalism because we are going to define the global expansion law as well as the 
large-scale smoothed trajectory field defined by averaging. In fact Buchert and Ehlcrs studied the averaging problem 
in Newtonian cosmology in the Lagrangian framework by performing spatial averages of Eulerian kinematical fields 
such as the rate of expansion 8 := V x ■ v of the fluid flow |I(],[l3|]. They have found that the background introduced 
by the spatial averaging obeys the FRW cosmology under the appropriate assumption that the peculiar velocity field, 
which is defined as a deviation from a Hubble flow in the Eulerian picture, obeys the periodic boundary condition on 
a sufficiently large scale (see below). 

In this paper we modify the approach by Buchert & Ehlers; we will work entirely within the Lagrangian framework. 
Namely, we divide the trajectory field into mean flow and deviation field, and then take the spatial average of the 
Lagrange-Newton system in order to introduce the horizon scale background as well as the large-scale averaged 
trajectory field. We arrive at the same conclusion as that of Buchert & Ehlers when averaged over the horizon scale. 
Then, in order to separate the non-linear dynamics from the large-scale dynamics, we further separate the deviation 
field into two parts, the averaged large-scale field and the rest. The evolution equation for the large-scale field is then 
obtained by averaging the local dynamical equation over a large domain much smaller than the present horizon scale 
in which the periodic boundary condition is applied for the small-scale perturbations. In this way we will obtain 
the evolution equations for the averaged large-scale field and the local-scale field. The evolution equation for the 
local-scale field naturally comes out by subtracting the averaged evolution equation for the large-scale field from the 
non-averaged equation. These equations couple to each other in general situations. Therefore, we are able to study 
how the smoothed large-scale structure is formed when the universe has non-linear structures on small scales as well 
as how the small-scale fluctuations grow in the surrounding environment. 

This paper is organized as follows. In Section 2, we shall write down the basic equations needed in our considerations 
in the Lagrangian formalism developed by Buchert. In Section 3 we investigate the average properties of the equations 
derived in Section 2. We will have the averaged FRW background under the periodic boundary condition over the 
horizon scale. In Section 4, we develop the formalism to have evolution equations for the large-scale and small-scale 
fluctuations where the results in Section 3 are used frequently. In general the evolution of the large-scale fluctuations 
is influenced by the existence of the small-scale nonlinearity. We clarify the situations where the large-scale fluctuation 
behaves independently of small-scale structures. In those cases the evolution of the large-scale fluctuations can be 
described by the so-called "truncated" or "optimized" Lagrangian perturbation theory which has been originally 
developed by many authors JlJ,[l7],^2],^3],[32| in order to avoid the shell-crossing problem of nonlinearity on small 
scales. Then it is also shown that the small-scale dynamics is governed by the modified scale factor. The final section 
contains discussions and summary. Throughout this paper, Latin indices take 1,2,3, respectively. 



2. BASIC EQUATIONS IN THE LAGRANGIAN PICTURE 



Let us start with the basic system of equations in Newtonian cosmology describing the motion of a self-gravitating 
pressureless fluid, so-called "dust" . The dynamics of the fluid obeys the following familiar Euler-Newton system of 
equations in Newtonian hydrodynamics: 



dp 



+ v x • O) = o, 



Of 

— + (v ■ V x )v = g, 
x g = 0, 

g = -4ttGp + Ac 2 , 



(2.1a) 

(2.1b) 

(2.1c) 
(2.1d) 



where p(x, t), v(x, t), a nd gj x, t) denote the fields of mass density, velocity, and gravitational acceleration, respectively. 
The Poisson equation ( 2. Id ) is extended including the cosmological constant for the sake of generality. 

Following the Lagrangian formulation developed by Buchert we concentrate on the integral curves x = f(X, t) 
of the velocity field v(x, t): 



df 

dt 



f 



df 

dt 



«(/,*), 



f{X,ti)=X, 



(2.2) 



A" 
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where X denote the Lagrangian coordinates which label fluid elements, x are the positions of these elements in 
Eulerian space at time f, and tj is the initial time when Lagrangian coordinates are defined. 

Then we can expr es s th e fields p, v and g in the Eulerian picture in terms of the Lagrangian coordinates (X,t) 
H| from Eqs.(^Ta|), and (|2.1b|) , respectively: 



p(X,t) 



Hx) 



(2.3a) 



j(x,ty 

v(X,t)=f(X,t), (2.3b) 
g(X,t) = f(X,t), (2.3c) 

where J is the determinant of the deformation field f^ij (the vertical slash in the subscript denotes partial derivative 
with respect to the Lagrangian coordinate X) and the quantities with ° such as P denote the quantities at the initial 
time ti henceforth, and we have used the fact J— 1. Thus, the continuity equation (2.1a) can be exactly integrated 
along the flow lines of the fluid elements in the Lagrangian picture rajqj9|. As a resu lt, the dynamical variable in the 
Lagrangian picture is only the trajectory field /. The equations ( |2.3b|) and (2.3c) are similar to point mechanics. 
The constraint equations ( |2.1c ) and (2. Id) of the acceleration field g give us the four evolution equations of the single 
dynamical field / after the usual procedure in the Lagrangian formalism: 



£abcfi\afj\bfj\c 0, 
n^jab^icdfc\afd\bfi\j 



Ac 2 J = -4vrG P (X), 



(2.4) 
(2.5) 



where we have used the mass conservation ( 2.3a ). The set of equations can be solved in principle for some initial 
conditions if we give the ini tial d ensi ty fi eld P( X) a s a source function. 

The system of equations ( ^.3a| -c), (2.4), and (2J5) is the so-called " Lagrange-Newton" system and equivalent to the 
Euler-Newton system as long as the mapping / f : X ^ s is invertible. Buchert has solved the above set of equations 
perturbatively taking the solutions of FRW models as the zero-th order background solution (as discussed below) 

m 

For the purpose of the later discussion and as an illustration of the Lagrangian formalism, we consider the Eulerian 
vorticity fiel d ujj = (l/2)eijkVk,j and derive the Kelvin's circulation theorem in the Lagrangian representation |l|,^,^|. 
First, from (2.3b) we can rewrite the vorticity field u> in terms of the trajectory field / as 



2,1 



£abcfj\afi\bfj\ 



To derive the theorem, we need Eq.(|2.4|) 
obtain 



This equation can be rewritten as 



Multiplying the equation by the hi 



£lbcfj\bfj\c — 0. 



dt \~"»"3\° J S\ 

so it can be integrated exactly along the trajectory field 

1 
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j t [eibcf 3 \bfj 

>ry field: 

£blcfj\bfj\c =UJ i, 



(2.6) 

the inverse matrix of f^ a , we can 

(2.7) 

(2.8) 
(2.9) 



where we hav e used the initial condition Wj= ( 1 /2 ) ej jk fh \j ( X , t j) from Eq.([2.6]). Finally, multiplying the above 
equation ( |2.9| ) by the deformation field fi\i/J and using Eq.(|2.6D, we can obtain the following Kelvin's circulation 
theorem along flow lines we are looking for: 



ulX,t) = jf ilj Z j (X). 



(2.10) 



The vorticity field in the Eulerian picture evolves according to the above equation along the flow lines and is coupled 
to the density enhancement, because th e field is proportional to the inverse of the determinant of the deformation 
field as the mass conservation equation ( 2.3a ) ||. This equation also means that if the initial vorticity field is zero 
at some point, the vorticity field remains zero at any later time along the flow lines. Conversely, as the density field 
develops singularities (J — > 0), the vorticity field will blow up simultaneously even if the initial vorticity field is much 
smaller than the irrotational part and is not zero [^). Therefore, one should bear in mind that the vorticity field might 
play an important role in structure formation in the non-linear regime. 
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3. AVERAGING NEWTONIAN COSMOLOGIES 



A. A Hubble flow for a trajectory field 



Before discussing the averaging problem, we first consider the properties of the trajectory field / in the Lagrange- 
Newton system derived in the previous section. We have reduced the description of the dynamics of any Eulerian 
field to the problem of finding the field of trajectories / as a solution of the Lagrange-Newton system (2.4) and 



(2.5). As in the Eulerian case, we are not able to write down any exact solution for generic initial data P (X) 
without assuming a symmetry like plane or spherical symmetry. We may start with the simplest class of solutions, 
the homogeneous-isotropic ones, and then move on to the treatment of inhomogeneities. 

Those fluid motions which are locally isotropic in the sense that, at any time and for each particle V, there exists 
a neighborhood on which the field of velocities relative to V is invariant under all rotations about V, are given by the 
following form with our choice (2.2) of Lagrangian coordinates: 



x = f H (X,t) = ^-X, S:=a(t/), 
a 



(3.1) 



if we conventionally put 
Lagrange-Newton system (2.4) 



2J), (U) ar 



Such a flow is a well-k nown Hubble flow. Inserting this ansatz into the 
and the mass conservation (2.3a) yield the usual Friedmann equations 



oft) 
a(t) ' 

Pb(t) 



4nG P b a 3 Ac 2 
" + IP 



3 a 3 {t) 



o3 

a 



Pb 



(3.2) 
(3.3) 



Thus, the quantity a(t) agrees with the scal e fa ctor in FRW cosmology. It should be remarked that the assumption 
of homogeneous and isotropic matter flow (3.1) makes the initial density independent of X via the equation ( |2.5| ): 
P{X) =Pb— constant. In o ther words, the existence of the fluctuation for the initial density field produces no longer 
a Hubble flow such as (3d). We may use it to integrate Eq^^) yielding Friedmann's differential equation: 



K.c 2 _ 8irG P b a 

-,2 ~~ 



o3 



Ac 2 



ft- 



ft' 



3 + 3 



K, = const. 



(3.4) 



where K, is the consta nt o f integration mathematically and can be regarded as the curvature parameter of the FRW 
model. Naturally, Eq.(2J3) gives 



u = 0. 



(3.5) 



Thus the assumption fl3.l|) for the trajectory field produces the standard Friedmann cosmologies. 



B. Averaged properties of the Lagrange- Newton system in an inhomogeneous universe 

We now consider the trajectory field when there exist inhomogeneities in the universe. In the application of 
Lagrangian theory to the averaging problem in cosmology, we examine the behavior of some spatially compact domain 
V(t) on the Eulerian space occupied by the fluid elements, which corresponds to the initial domain T> of the Lagrangian 
coordinates via the mapping f t : X i— » x. 

For our purpose, we set the average flow in the form of a Hubble flow with scale factor o-p not necessarily equal to 
a(t) and define the (not necessarily small) deviation field P from the average flow without loss of generality, so that 
the full trajectory field / of an inhomogeneous model reads: 

f(X,t):=f^X,t) + P(X,t) = ^-X+P(X,t); P(X,tj):=0, (3.6) 

Thus, for the sake of convenience, we here start with the trajectory field in the form of Zel'dovich type solution 
Q. However, it should be noted that we can define the full trajectory field by Eq.( |3.6|) for any fluid elements in an 
arbitrary initial domain T>, and the scale factor depends naturally on the chosen domain. We also remark that since 
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we have imposed no condition on the deviation held P, this consideration is applied to the non-linear situations when 
the density contrast held may be lar ger than unity. 

Using the full trajectory held (3.6), we obtain by a little computation the following expression for the determinant 
J of deformation held f^u: 

J 7y£abc£ijkfa\ifb\jfc\k 

= 7T + 7T P *\i ( P i P ]\j ~ P i P i\j)u + 7 [ P i { P j\3 P k\k - p j\k p k\j) + 2Pj {Pi\k p k\j - P k\k P i\j)] u 

a v a v z a T> 1 o 

:=±-(aUt)+Ji\i(X,t)), (3.7) 



where 



Ji -=av alPi + ^-a v (P^ - P,P tb ) + ^ [P t (P 3 yP k \ k - Pj\ k Pk\i) + 1 p i { P i\k p k\j - p k\k p i\j)] ■ (3-8) 



3 2 

27 _ (r, n n n ^ , "27 



Note that the second term on the right-hand-side of Eq.(3.7) is expressed by the divergence of the vector jT(X,t) 
with respect to the Lagrangian coordinates and it is defined in terms of the deviation displacement vector P. Since 
the volume elements at t and ti are related by d 3 x = Jd 3 X , we use the above equation to rewrite the volume Vvif) 
of the domain T>(t) of the fluid in the following form: 



V v (t):= [ d 3 x = fd 3 XJ(X,t) 
Jv(t) Jv 



1 



a 



v 



Ji\i) )Vv 



a \vv+\ I dS x -J(X,t), (3.9) 



where we have applied Gaufi's theorem to transform the volume integral to a surface integral over the boundary dT> 
of the initial domain: 

r d 3 XW x J= fjSx J- 

27 JdT> 



In the derivation of Eq.(3.E), (. . .)^ denotes the spatial average of a tensor held over the initial domain f), and we 

regard a-p(t) as the scale factor of that domain V(t) in the Eulerian space. The quantity Vv denotes the volume of 
initial domain T> considered, 



VV-= / d 3 X. (3.10) 



^3- 

'27 

A note of caution is in order: Buchert & Ehlers (h]] have used a domain dependent scale factor av defined by 
Vv = dj){t), but Eq.( |3.9| ) means that such a scale factor does not agree with our scale factor dehned by Eq.( [j.6| ). 
Thus, these two concepts arc different in the following respect: if we would impose periodic boundary conditions for 
P on the domain, then our scale factor reduces to t he sta ndard FRW scale factor on that domain as discussed below. 

Likewise, the spatial average of the density held ( [2.3a ) over the domain V in Eulerian space at any time t may be 
calculated as follows. 

i i J3 . !m 



d 3 Xp 



1 

Vv 
Vv 

Vv(t) V '27 

''V 



■(Ph 



o3 



dSx-J 

dt) 



{R{X)) h . (3.11) 
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Next we consider t he a verage of the dynamical equation (£^j) of the trajectory field. Likewise, the first term on the 
left-hand-side of Eq. ( |2.5| ) can be rewritten by inserting Eq. ( |3.6[) in the following form: 



where 



7: e abc£ijkfa\ifj\bfk\c — — (3d|,ap + Q i | i (X,t)) 

1 a-v 



(3.12) 



+ a v a v [P t P jU - P 3 P A3 ) + -fd v (PiP^ - PjPy) 



o3 



P i \ P j\j P k\k - P]\kPk\j) 



2Pj (Pi\kPk\j — Pi\jP> 



k\k, 



(3.13) 



and we again remark that the second term on the right-hand-side of Eq.( 3.12 ) can also be expressed by the divergence 
of the vector Q(X, t). Thus, we can get the local evolution equation for the domain dependent scale factor a-p and 
the displacement vector Pf. 



.3 



Sa^a-D + 47rG P{X) a v -Ac 2 a^ 



i | i 



Ac 2 J, 



(3.14) 



Averaging over the initial domain T> of the above equation (3.14) leads to the following equation: 



3^ + 4ttG 



' v -Ac 2 



dS 



x 



Q + Ac 2 J) , 



(3.15) 



where we again used Gaufi's theorem. If the local equation (3.14) and the averaged equation ( |3.15 ) can be solved 
simultaneously, the domain dependent scale factor az>(t) and the local displacement vector Pi are obtained in principle, 
respectively. This averaged equation can also be interpreted as a standard Friedmann equation for the "effective mass 
density" p c g IhJ , which is here defined by 



4irGp eS (t) 



i ° , o 3 

4irG — ^ 



dS 



X 



Ac 2 J 



(3.16) 



where the first term on the right-hand-side decreases clearly in proportion to a-p 3 

inhomogeneities have an accelerating effect on the expansion rate dp /a-p of the average flow, if the term J^a dSx 



The equation ( 3.1 5| ) shows that 

(-C) 

on the right-hand-side dominates the other terms and is positive. Namely, this shows that the evolution of the domain 
dependent scale factor a-p does in fact depend on the chosen domain, that is, the averaged expansion will be different 
from the usual Friedmann laws (ET^) if the averages involving Q and S do not vanish. 



Next, we consider the property of the average of the vorticity field u>. Similarly, performing the average of Eq. ( 2. 1C ) 
over the domain T> on the Eulerian space (not the Lagrangian space), we can obtain 



o3 



i 3 v (t) + (1/ 



Vv / dSx ■ J(X, t) 



up 



a o 



1 



Vv J 9T> 



dS 



Pi u. 



(3.17) 



Noting that the quantity w.j is divergence-less by definition, namely ^1^= (1/2)6^^ v^.i— 0, we have again used 

Gaufi's theorem in rewriting the second term on the right-hand-side. The equation (3.17) means that if cij vanishes 
at every point on the initial hypersurface, it leads to (lo^j, = at any time. 

Thus, the average properties of the Lagrange-Newton system do not necessarily agree with the FRW cosmologies 
in a general inhomogeneous universe. However, based on the observation of extreme isotropy of the CMB, we expect 
that the universe is almost isotropic and homogeneous on a sufficiently large scale. We consider how this fact is 
expressed mathematically in terms of averaged variables. 

As discussed above, we could define the displacement vect or P as representing the deviation from the mean flow 
generated by the inhomogeneities. The resultant equation ( 3 . 1 5| ) then shows us how this field P determines the 
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backreaction on the scale a-p(t) from Friedmann's law. Note that the backreaction terms are expressed by the surface 
integrals over the boundary of the initial domain &D- As already noted, if we employ the periodic bo unda ry condition 
for the deviation vector field P on some sufficiently large scale T> p , the backreaction terms in Eq.( |3.15 ) are exactly 



AttG 



o o 3 

p b a v 



Ac 2 = 0, 



(3.18) 



with 



p b = (p (-X"))o, = constant. 



Similarly, Eq.( 3.11 ) gives the background density at an arbitrary time: 



Pb (t) = (p(X,t)) 



O o 3 

p b a v 



v p (t) 



.(*)" 



(3.19) 



(3.20) 



We thus obtain the usual definitions of the homogeneous background density field ( 3.19] ) and ( 3.20| ) using the spatial 
averaging. Since we have not restricted ourselves to the perturbative situation where the deviation vector P is 
infinitesimally small, this discussion is always valid for non-linear situations under the periodic boundary condition. 

Under the same periodic assumption, Eq.(3.17) gives us the following form as the average of the Eulerian vorticity 
field: 



(3.21) 



with the scale factor ar> p defined by Eq.(3.18). Thus the averaged vorticity field decays as oc a v 2 in the expanding 
universe, which is analogous to the Newtonian linearized theory |p6| . In other words, even if the initial global averaged 
value (wi)^ is not zero, we can safely ignore the global averaged vorticity field. 

The COBE microwave background measurement suggests that the power spectrum of the density fluctuation field 
has a positive slope on large scales, supporting the assumption of large-scale homogene ity. This suggests that even 
if we do not employ the periodicity over the horizon scale, the flux of Q and in Eq.(3.15) through the boundary 
of the averaging domain with a sufficiently large volume may be negligible. Therefore, we may conclude that the 
backreaction on the global expansion rate becomes zero and the equations of the background mo del introd uced by 
the spatial average of an inhomogeneous universe over the horizon scale obey Friedmann's laws (|3.2|) and (3^3) in 
Newtonian cosmology, even when the universe ha s lo cally nonlin ear structures (5 ^> 1) on some small scales. 

We may proceed to solve the set of equations (^4) and ( 3.14 ) perturbatively by taking the solutions of the Fried- 
mann's laws as the zeroth order approximation for the scale factor to construct a locally inhomogeneous universe 
[fifrl . This approach is certainly useful, but it is not easy to see the effect of large-scale structure on the small-scale 
nonlinear dynamics. To see this explicitly we will take a new approach in the next section. 



4. HYBRID LAGRANGIAN THEORY FOR NEWTONIAN GRAVITATIONAL INSTABILITY 

In 1970 Zel'dovich |3^] found the so-called Zel'dovich approximations to describe the large-scale structure formation; 
it has been derived in the gravitational context by Buchert Q and it can also be obtained as a particular first-order 
solution in the Lagrangian perturbation theory developed by Buchert and many other authors [^, p^ , p9[ . It 

has been shown that Zel'dovich approximation gives indeed a very good approximation and has been used up to the 
quasi non-linear regime to reproduce the observed filament-like and pancake-like pattern of the large-scale structure 
beyond several megaparsecs. 

The numerical simulation based on the approximation has an advantage over the N-body simulation such that it 
is able to simulate relatively large domains with relatively small memories. However, the Zel'dovich approximation 
also has a disadvantage. Namely, it cannot reproduce the non-linear structure formation on small scales once the 
shell crossing occurs. This is called the shell crossing problem. To avoid this, the truncated or optimized Lagrangian 
perturbation approaches have been developed by many authors jl2],[l7 22 2^j3^] where one smoothes out the initial 



7 



small-scale fluctuations such that the shell-crossings occur at about the present time. This will be a good approxi- 
mation for the evolution of the large-scale structure, if one can show that the evolution does not very much depend 
upon the behavior of small-scale nonlinear dynamics. However, this expectation has not yet been explicitly proven 
or, put in another way, it has not yet been clarified under what sort of situations this expectation is valid. If this is 
proven, there may be some hope to have a hybrid way to describe the inhomogeneous universe such that the large-scale 
structure is described by Zel'dovich's approximation or some improved version of it and the local small-scale dynamics 
is described by other method such as N-body simulation or some effective theory of nonlinear structure formation. 

This is what we wish to develop in this section. We do this by dividing the deviation field into two parts, the 
averaged large-scale part and the small-scale part. In the following, we consider the Einstein-de Sitter universe 
(A = K, = 0) for the sake of simplicity, and the scale factor is normalized as 

a B = 1. (4.1) 
Namely, we set the scale factor so that the comoving coordinates agree with the Lagrangian coordinates. 



A. The division of the deviation field into the large-scale part and the small-scale part 

Let us start writing the trajectory field as 

MX, t) = /f (X, t) + Pi{X, t; A) - a(t)Xi + P^X, t; A), (4.2) 

where we have explicitly included the wavelength (A) dependence in the deviation field Pj. Note that A here denotes 
the wavelength of the initial density fluctuation field in the comoving coordinates as explained later. Here the quantity 
a(t) is the scale factor defined in the previous section: a = ap , (where we have omitted the subscript of the quantity 
a-p for simplicity). Namely, the scale factor is obtained by averaging of the Lagrange-Newton system over the 
horizon scale. Following the results in the previous section, we assume henceforth that it obeys Friedmann's law, so 
the deviation field Pi obeys the periodic boundary condition on the present horizon scale: 

(P A] {X,t-K\<L H )) a = 0. (4.3) 

v > J i Vh 

The wavelength of the initial density fluctuations has the lower cutoff, because we deal with a collisionless gravitational 
system like dark matter. Below the lower cutoff this description is not valid since the baryonic gaseous pressure and 
the effective pressure due to the velocity dispersion of the collisionless system may become important. We will not 
consider such small scales. We have also the upper cutoff which will be the horizon scale Lh because we assumed that 



the deviation field Pi obeys the periodic boundary condition (4.3) on the horizon scale Lh- We will not explicitly 
write down these cutoff lengths hereafter. 



By substituting the ansatz (4.2) into Eq.(2.5), we can obtain Eq.(3.14), or more explicitly as 



3a 2 a + a 2 Pm + 2aaP,;| i + a (Pj , / ', , — P i\j P 



+ ^ ( p i\i p j\j ' p i\j p j\i) + 2 e ^ e abcPi\ a p j\b p k\c = -4ttG P {X- 1 < A < L H ). (4.4) 

Now we introduce the averaged vector field pf by using the spatial average of the full deviation field Pi over the 

large-scale domain T>l(X) (which is still much smaller than the horizon scale) at some point X in the Lagrangian 
coordinates: 



p>(X,t;\>L) = J- ( (PX'PiiX'^X) 
Vl Jv l {X) 

:= -!- / <PX'P t (X', t; X)W(X - X'; L), (4.5) 
Vt. J 



where L is an artificial cutoff length and 



Vl:= L d 3 X= [ d 3 X'W(X-X';L), 
Jv L {x) J 



(4.6) 



and W(X — X'; L) is a filter function characterized by the smoothing length L. For example, for Gaussian filtering 
on scale L the filter function is 
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W G (r;L) = 



( 27r )3/ 2jL 3 

From Eq. fl3.6[ ) the initial condition is imposed on the pf : 

P >(X,t!) = 



exp[-r 2 /(2£ 2 )]. 



(4.7) 



(4.8) 



We have assumed that, since the smoothing method with the scale L is used for the definition of pf , its valid 
wavelength range has about the length L as a lower cutoff length. Thus, we expect that the cutoff scale length L 
is large enough so that the pf describes the evolution of fluctuations with the characteristic length larger than L in 
the linear regime still at the present time to, that is, \pf^(X, to) | -C a(to) is satisfied according to the Lagrangian 
perturbation theory, ft should be noted that the linear regime in the Lagrangian picture does not mean j< 1 because 
the density field can be exactly solved. In this situation, we can safely say that the cutoff length L is supposed to be 
in the range 1 « I « Lh in the realistic universe. In the following, we consider only the lowest order in pf by means 
of this assumption. 

Simultaneously, since the original displacement vector Pj d escribes the gravitational instability of fluctuations with 
scale larger than I and smaller than Lh, the definition (L5) of pf allows us to define the small-scale displacement 
vector pf which is supposed to describe non-linear structure formation on small scales. Namely, the full displacement 
vector Pi can be written as 



Pi(X, t; A) = pf (X, t;L< 1 X) + pf (X, t; A £ L). 



(4.9) 



The initial condition for pf is 



P f(X,t I )=0. 



(4.10) 



We expect that pf describes the evolution of the initial density fluctuation field with much smaller length than our 
cutoff length L. 

Next, we shall derive the evolution equations for the large-scale displacement vector pf and the small-scale dis- 
placement vector pf , res pectively. Before doing this we present a basic equation which is obtained by inserting the 
ansatz (lOjj) into Eq.(O). 



3a 2 d 



^P% 



+ ' d W\iPf\i 



^p% 

pf .p.,. 



2aapf. + 2ap" 



(p?\iP 



< 



i\i 



pfPk\ k -PkPf\k 



P 3 Pk\k -PkPj\ k 

\pfPk\X k -(Pk\kP. 



Pi\jPj\i 
2 P%Pt\. 



(p%p 



Pk | kPj | i 



< 



P%Pf\i 



+a 



Pfp% 



pfp K 



>\2\ - 



L ° ( < < 

2 \ Pi P M 

+0({ P > ) 2 ) =-^Gp{X;l<X<Lu). 
In the above derivation, for example, we have used the results such as 



Pk\kP% 



Pf\kPk\j Pf\jPk\k 



(4.11) 



Pi\iPj\j 



Pi\jPj\i 



PiPf\j-Pfpf\s 



Note that we have kept only linear order in pf and full order in pf . It should be also noted that only the third term 
in the second line on the left-hand-side of Eq.(4.11) cannot be expressed in the form of a divergence of the vector 
which consists of pf. In the discussion below, the large-scale transverse mode is omitted in the sense that we use 
only one of the equations (2.4) and (2.5), namely (4.11). For the smoothed large-scale field pf, this may be a good 



approximation, because assuming that the initial vorticity field is negligible compared with the initial irrotational 
flow results in vanishing of the transverse part at a later time as explained in Appendix A (this may be a reasonable 
assumption based on the linearized theory [2f|), but there are transverse parts in the Lagr angia n space in the nonlinear 
situation even for uj 1 = 0. Buchert & Ehlers || derived the transverse solutions of Eqs.( [2.1a| -d) for the second-order 
transverse and irrotational solutions in a general case. 
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B. The evolution equation for the large-scale structure formation 



First, we consider the evolution equation for the averaged large-scale field pf . Since we have assumed that the pf 
describe the behavior of fluctuations with the characteristic scale much smaller than L, we introduce the following 
spatial averaging method on the scale V much smaller than L and much larger than I: 1 < L' < L. Namely, the 
averaging is defined by 

(1) (g) L , (X, i) s -L f d s X'g(X',t) := J- / d 3 X'g(Xf t)W(X - X'; L>). (4.12) 

Vl> Jv l ,(x) Vl> J 

Furthermore, within the averaging volume V v-, we can safely neglect the spatial gradient of pf because from Eq.( fl.5| ) 

the pf are defined by the averaging on the volume Vl with L 3> L' . Thus, we can safely employ the following second 
rule when we perform the averaging of an arbitrary function F (p < (X , t) , t) of the small-scale fluctuation field p < 

multiplied by p > over the volume Vl 1 '- 

(2) -J- f d 3 X'pf(X',t;L<\)F(p<(X' 7 t),t)=pf(X,t;L<\)^- f d 3 X / F(p<(X' ,t),t). (4.13) 

Vl> j v l , (x} Vl< Jv l ,(x) 

This rule is correct at the lowest order in a Taylor expansion of p > : 

pf{X\t) = pf(X, t)+p> j (X, t) (X - X). + ■ ■ ■ , (4.14) 

because we have neglected the second term in the above compared with the first term. In the averaging volume Vl', 
the first and second terms are of order p > and (p > L')/L, respectively, so the assumption L' <C L allows us to ignore 
the second term compared with the first term. Finally, we introduce the third rule for an arbitrary vector Gi(p < ) 
which consists of the small-scale displacement vector pf : 

(3) 4- I cPx'Gi^) := [ d 3 S' l G l (p<)=0. (4.15) 

Vl> Jv l ,(x) Vl< JaD h ,{x) 

Here, we have assumed that the p < are mainly generated by random initial density fluctuations with only the char- 
actcristic wavelength much smaller than L , and obeys the periodic boundary condition on the volume Vu v^-V l)'- 



4i)r, = ^~l d 3 X'p<(X',t)=0. (4.16) 

U/ L' Vl>Jv l ,(X) U 

This third rule causes a possible error because we neglect the fluctuations with scale L' < A < L. Neglecting the 
fluctuations with the length comparable with the scale L may not be a serious problem, because we have assumed that 
the fluctuations with the scale L are still in the linear regime at the present time. However, neglecting the fluctuations 
with scales comparable with L' might cause a serious problem. For the present we leave this problem open, and we 
consider the situation under the above three rules. 

According to these rules, by averaging both sides of Eq.(l.ll) over the domain Vl 1 we obtain 

3a 2 a(t) + a 2 p%(X,t) + 2aap>. + P > ]j (X,t) (p<.p< k - p<.p< k ) ^X,t) = -4nG(°p) L '(X- L' £ A), (4.17) 

where we have used the following calculation 

( 2 Pf\kPk\i ~Pf\iPk\k -pf\iPk\k) v = I 2 (pfPk\i) ]k - 2 Pfp%k ~ Pf\iPk\k ~ Pf\iPk\kj t 

2 (pfPk\k) {i + 2 Pf\iPk\k - Pf\iPk\k ~ Pf\iPk\k ; 



(pf\iPt\k-Pf\iPk\k) L ,( X ^- ( 4 -!8) 



Using the Friedmann equation ( |3.2| ), Eq.( 4.17 ) becomes 

a 2 p%(X,t) + 2aap>fX,t) +p> ]j (X ) t) (p< p< |fc - p< t p< k ) (X ,t) = -4nG ((P)l'(X;L> & A)- P b ) ■ (4.19) 
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We emphasize that the source function on the right-hand-side of the above equation became the density fluctuation 
field with a wavelength larger than the averaging scale L' and sm aller than Lh because of using the smoothing 



method. Furthermore, the third term on the left-hand-side of Eq.(4.19) represents a backreaction effect that the 
small-scale non-linear displacement vector p < has on the large-scale perturbation p > . Namely, even if the small non- 
linear displacement vector pf obeys the periodic boundary condition on the volume V l 1 , the non-linear structures 
have a possibility to give such a backreaction effect on the evolution of fluctuations with larger scales in the linear 
regime of the larger fluctuations. 

However, we have the following situations where the large-scale backreaction becomes zero or negligible compared 
with the other terms. 



(1) The first case is that the small-scale displacement vector pf can be divided into a time function and a spatial 
function with respect to the Lagrangian coordinates: pf — D(t)ipf. Then the backreaction term becomes zero: 



Pf\iPk\k ~P%Pk\k) L , = D(t) (^f\^k\k ~ ^f\i^k\k) 



= 0. 

During the early stage after the decoupling time of matter and radiation, we expect that the pf can be described 
by the Zel'dovich type solution and thus we may have the above case. Actually, Ehlers & Buchert Jl3| have shown 
that the displacement vector at every order has such a separable solution in the Lagrangian perturbation theory. 
But in the non-linear regime later around the peak patch of density fluctuation field, there is no guarantee that 
pf is separable. In this paper, we are interested in the non-linear situation of pf. 

(2) The second case is as follows. We can divide the deformation field of the large-scale displacement vector pf^ 
into the divergence, trace-free symmetric, and antisymmetric parts without loss of generality: 

ff\j = \ 5 aPk\k + Q(Piu +Pf\i) ~ l s »Pk\k) + \ipfy - Pf\i) = M a + <% + ^ ( 4 - 20 ) 



where 



A(X,t) = ±p>., (4.21) 



S l3 (X, t) = ~(p> + p> ) - -5 ijP > k , (4.22) 
MX,*) = i(p>-p>) 



2 X 

K l3 (X,t)^l(p>-p>). (4.23) 



As discussed in appendix A, we can safely ignore the antisymmetric deformation field IZij compared with the 
other quantities A and <Sy under an appropriate assumption: IZij w or 72-jj <C A, Sij . Then, if the expansion 
deformation A is assumed to be much larger than the shear part Sij , the backreaction term becomes 

P% {p%Pt\k ~ Pf\iPk\k) L , = My (p%p<\k ~ Pf\iPt\k) L , + S v {p%Pt\k ~ P%Pk\k) L , ^ O {(Sp<p<)) . (4.24) 



Hence, to ignore the backreaction term compared with the other terms in Eq.(4.19), we have to employ the 
following condition 

a 2 A, adA » Sy (pf^ k ~ P%pf lk ) £< ■ (4.25) 

This assumption must be handled with caution. Since p > represents the averaged large-scale field of the original 
displacement vector Pi, we expect that it can be expressed in the form of Zel'dovich type solution as discussed 
below. Thus, the assumption A ^> Sij may not be appropriate in the quasi non-linear regime, and we need 
more investigations in detail. The condition ( [4.25 ) could be checked in a realistic structure formation scenario 
by using the numerical simulation. 

• (3)The third case is that a locally one-dimensional motion or a spherical top-hat motion dominates for pf. in 
the nonlinear structure formation on small scales: 

Pf\j »P^i*«.*ii> or Pf\j™Pk\k 6 ij- ( 4 - 26 ) 
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In these cases, the backreaction term becomes 



pty {pfA - p%Pk\k) L , « E p% \pf\iPk\k - pf\iPk\k, v 

clump 

= J2 pf\i(pf\i\pf\i-pf\ipf\i) L , ° r J2 pt\i 5 ^{p%pt\k~pf\ l pt\k l 

clump clump 
= 0, (4.27) 

where Sclump denotes the sum of the number of non-linear small-scale clumps included into the averaging 
domain T>L' ■ Here we have assumed that we can replace the integral in (• • ■) L , with the sum over the clumps. 
• (4) Finally, we remark that the backreaction term can be rewritten as 

Pi\3 {p%Pt\k - Pf\iPk\k) L , = P%J t (pf\iPt\k ~ P%Pt\k) L> ■ ( 4 - 2 8) 

This suggests that we may employ a time averaging together with the spatial average to obtain the evolution 
equation of the large-scale dynamics. Once the nonlinear structures are developed, its time scale will be much 
shorter than that of large-scale dynamics. Thus, the averaging over the local time scale will eliminate the above 
backreaction term. Although this is interesting, we will not pursue it here and leave it for future study. 



Under these situations or some combined situation of them, the evolution equation (4.19) for the large-scale dis- 
placement vector pf yields 



a 2 p> + 2aap> + O ((p>) 2 ) = -4ttG P b 5 L '(X; L' < A), (4.29) 



where 



(p) L ,(X;L' <\)- p b 

6l>(X; U < A) = WLl - g ' -. (4.30) 

Pb 

We again note that, since pf is defined by smoothing the original displacement vector P,, we need only the lowest 
order of p > . 

Eq.(4.29) entirely agrees with the first order evolution equation in the Lagrangian perturbation theory, so we can 
proceed to solve it for pf iteratively using the solution of the scale factor a(t) for the Friedmann background equation 



( p.2|) . In this way, introducing the artificial cutoff length L and the spatial averaging over the volume Vl' smaller 
than the horizon scale allow us to understand the validity of the Lagrangian perturbation theory in investigating the 
large-scale structure formation even if the universe has non-linear structures on small scales. Our formalism clarifies 
the meaning of the optimized or truncated Lagrangian perturbation theory which has been frequently used. Many 
authors Jl2] , ^7| , p2y23|j32] ] have shown that the optimized or truncated Lagrangian perturbation theory reproduces the 
results of the large-scale structure larger than the smoothing scale by the full N-body simulation. This indicates that 
some of the above assumptions may be satisfied in the situation. It may be interesting to investigate which assumption 
is valid, and this will be presented elsewhere. 



C. The local small-scale non-linear evolution equations in the Lagrangian picture 

Next we derive the local evolution equation for the small-scale displacement vector pf . We have obtained the back- 
ground equation and the evolution equation for the large-scale displacement vector within the Lagrangian framework. 
We want to construct the local evolution equation for the small-scale displacement vector including the effect of the 
gravitational instability of the surrounding large-scale structure. If we do so, we may have the possibility to use 
the high-resolution N-body simulation or the semi-analytic approaches only for the small-scale non-linear structure 
formation and to use Zel'dovich's approximation for the large-scale structure formation, simultaneously. 

Subtracting Eq.(4.19) from Eq.(Lll), we obtain the following local small-scale evolution equation in the situations 
taken up in the previous subsection: 
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a 2 V% + 2a«4 + a (p> p<. - p>.p<^ + a - p> p<<) + a (,,] 



{p?\iP. 

( 



< 



Pf\jPf\i) + p% (I'i /'x I'l ji'ji ) +Pi\i 



< P < r 

> 



0. 



■p%p% 



j\kPk\i 



2 



< 



Pi\jPj\i 



' Pk\kPj\i 



+ 9.P% [ P f\jPk\k ~ Pk\jPf\k) +P% (Pf *■ 



I < < - 



fcPfcji Pk\kPf\i 



-4irG(°P) L , 5(X;\&L'), 



where 



P(X;A)-(P) £ >(X;L^A) 
(p) L ,(X;i'< A) 



(4.31) 



(4.32) 



Note that the source function is a density fluctuation field with scale smaller than L'. 



If we could solve Eq.( 4.29| ) for the pf , we could solve Eq.(4.31) for the local displacement vecto r pf in principle 



by substituting the solution p > . However, in practice, it will be very difficult to do so because Eq.( 4.31 ) is a highly 
nonlinear differential equation. Instead we restrict ourself to a more simple situation in this paper in order to see 
clearly the environmental effect on the small-scale dynamics. Namely, we consider the second situation described in 
the previous subsection. In this case the Lagrangian divergence of the displacement field dominates on the averaged 
large-scale dynamics. It is straightforward to extend the formulation taking into account the effect of tracefree part 
of the large-scale deformation field which represents the surroundi ng ti dal field, and it will be presented elsewhere. 
Here, we consider only the simplest case. Then we can rewrite Eq.(4.31) as 



(a + A) 2 p< {l + 2(a + A) (3 + A)p< +(a + A) (. 



PfuPf 



i\i F 3\i 
1 

Or, using Eqs.( 3.1§l ), ( 4.29 ) and ( 4.30| ), this equation becomes 



Pf\jP% 



+-(a + A) U ,//, , - ,/' ; , J + ^ ljk e abc p< a p< b p< c 



AnG(°P} L , °6(X;\&L'). 



(4.33) 



3 (a 



PfuPf 



A) 2 (a + A) + (a + A) 2 pf H + 2(a + A) (a + A)pf H + {a + A) (/*, j>\ , - 

+-(a + A) (p%P% -P?\jPf\i) + 2 e v^Pf\aPf\bPk\c = "4ttG P(X; A), 



(4.34) 



Note that we have normalized the global scale factor to a= 1 in this secti on a nd considered the Einstein-de Sitter 
background. If the above equation is compared with the original equation (4.4) in the Lagrangian picture, it is seen 
that the effect of the large-scale structure formation on the small non-linear scales is represented as a modification of 
the global scale factor (a + A)(t) on the surrounding large scale. 

Thus, since the spatial gradient of the expansion A with respect to the Lagrangian coordinates can be ignored 
inside the volume V v m this case, the above equation corresponds to the following system of Eulerian equations on 
the comoving coordinates y = x / (a + A): 



dp n a + A d ,. 

-dt +i ^TA p + W {pu) = ^ 

du l a + i . 7 - du* _ 1 

dt + a + A U+U dyi~ (a + A) 2 W l1 
A y( p = AnG{a + A) 2 ( P (y, t) - (p) L ,(t)) , 



(4.35a) 

(4.35b) 
(4.35c) 

where u 1 denotes the peculiar velocity on the comoving coordinates {y} and the physical peculiar velocity is (a + A)u l , 

and (p) L i(t) represents the effective background density defined by (p) L ,(t) = (p) l> / {a+A} 3 . The comoving coordinates 
{y} should be defined on the Eulerian coordinates at the initial ti me, wh en fluctuations with all scales are much 
smaller than unity. Thus, we can interpret the system of equations (4.35a-c) as describing the evolutions of small- 
scale fluctuations on the effective FRW background model characterized by the modified scale factor (a(t)+A(t)). This 
equation means that, because of the fact A < within the large-scale objects that are collapsing, the fluctuations on 
smaller scales tend to collapse earlier than the fluctuations in the large-scale void where one has A > 0. Furthermore, 
we can perform the N-body simulation using the above Eulerian set of equations inside the box with volume V v in the 
usual way only if the scale factor is modified to the effective scale factor (a + A) inside the box by using the Zel'dovich 
solution or the improved version of it for A(t). In this way, our formalism will be useful for the consideration of the 
environmental effects on the behavior of the fluctuations. 
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5. DISCUSSION 



We have developed a formalism which allows us to investigate the relation between large-scale quasi-linear dynamics 
and small-scale nonlinear dynamics using the averaging method in the Lagrangian theory. We have derived the 
coupled equation for the large-scale dynamics and the small-scale dynamics. In the case where the averaged large- 
scale dynamics is expansion dominated, we have shown that the large-scale dynamics decouples from the small-scale 
nonlinear dynamics. Then, on the other hand, the small-scale dynamics is influenced by the large-scale dynamics in 
such a way that the local small-scale equations contain the modified scale factor of the large scale. The modified scale 
factor is the sum of the global scale factor and the expansion of the region considered. Our result strongly suggests 
that there will be more complicated environmental effects in local small-scale dynamics which one cannot ignore. 

There may be several possibilities to generalize our analysis. One would be to employ some approximation to solve 
the local dynamics in the Lagrangian framework. For example the spherical symmetric approximation on small scales 
may be reasonable and we may employ the Press-Schechter |27| type approach to the local region |^l|,|^,^5). More 
challenging would be to employ N-body simulation on small scales. We have mentioned this possibility in the previous 
section, but it seems more work is necessary to achieve this consistently within our scheme. One thing missing in our 
formalism is the mutual gravitational interaction between the small clumps contained in a large scale environment. 
This is because we have ignored the dynamical freedom bet ween L 1 < A < L. Also, although we have here ignored 
the large-scale tidal effect on the small-scale dynamics in Eq. ( 4.35a -c) for the purpose of a simple illustration, it may 
play an important role in the hierarchical structure formation. For example, Bond & Myers Q have investigated 
the important influence of the tidal effect for the merging history of halo objects. It is straightforward to extend 
our formalism including the large-scale tidal effect on the small-scale dynamics, and investigating the environmental 
influence will be interesting. We would like to come back these points in future study. 
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APPENDIX A: THE LARGE-SCALE VORTICITY MODE IN THE LAGRANGIAN PICTURE 



In this appendix, we investigate the behavior of the rotational deformation field of the large-scale displacement 
vector pf . For this purpose, it is convenient to make use of the Kelvin's circulation transport equation (2.9) in the 
Lagrangian picture: 



2 €blc f]\bfj\c —^l ■ 

By substituting the full trajectory field fi — aXi + pf + pf into the above equation, we obtain 
1 



*(Pii* + Pj\k) - "(Pj\k + Pf\k) + Pl\jPf\k + Pi\kPi\j + Pi\jPf\k — ~ w i {X; I < \ < L H ) 



(Al) 



(A2) 



If we perform the averaging of the above equation over the domain Dv according to the rules ( 4.12| ), (4.13) and 
(|4.15p, we get 



I d Pj\k 
2 eijk dt a 



(A3) 



In the above derivation, for example, we have used the calculation such as 

-eijk {p{\jPf\k) L , = \^jk ((P(\ 3 Pf )\k) L , = 0. 



(A4) 



By noting pf^(X,ti) = and using the variable IZij defined by Eq.(4.23), Eq.(A3) can be integrated as 



-a(tl>i)i 



dt 



u a\t>) 



(A5) 
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Thus, if we assume that the averaged initial large-scale vorticity field {u>i) L , is exactly zero, we arrive at the conclusion 

Kij = 0. (A6) 

According to this consideration, even if we do not adopt the above assumption, we can safely ignore the large-scale 
rotational field IZij compared with the trace part of the deformation field A and the trace- free symmetric part Sij, 
because the initial vorticity field is much smaller than the expansion field and the shear field based on the linearized 
theory ||. 
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